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We suggest a new scenario in which the Universe starts its evolution with a fractal topological 
structure. This structure is described by a gas of wormholes. It is shown that the polarization of 
such a gas in external fields possesses a spatial dispersion, which results in a modification of the 
Newton’s law. The dependence on scales is determined by the distribution of distances between 
throat ends. The observed in galaxies logarithmic correction confirms that the distribution has 
fractal properties. We also discus the possibility of restoring such a distribution from observations. 


I. FRACTAL TOPOLOGICAL STRUCTURE AS 
INITIAL CONDITIONS 

The most important problem of the modern astro¬ 
physics is the nature of dark matter (DM). While particle 
physics suggests too many speculative candidates, LHC 
does not provide even any hint for the existence of such 
dark particles. Moreover, the standard model in particle 
physics needs not their existence at all. This make us 
think that something wrong with our interpretation of 
the picture of the Universe observed. 

It turns out that lattice quantum gravity, has already 
suggested a satisfactory principle answer to the question 
what DM should be. The answer lies in the very com¬ 
plex spacetime foam picture which exists at very small 
(sub-Planckian) scales and which possesses fractal prop¬ 
erties. Indeed, the first indication that at very small 
scales spacetime has a universal fractal structure came 
from 2d-quantum gravity [l]. This result was obtained 
both, analytically and by means of numerical simulations 
and it has been developed in many papers, e.g., see 0 
and references therein. Moreover, numerical simulations 
demonstrate that the fractal behavior retains in 4d lat¬ 
tice quantum gravity as well. And indeed, the spectral 
dimension for nonperturbative quantum gravity defined 
via Euclidean dynamical triangulations was calculated in 
Q. It turns out that it runs from a value of D= 3/2 at 
short distance to D = 4 at large distance scales. 

The fractal structure of the spacetime foam can be de¬ 
scribed as follows. Consider any point x in space, fix a 
geodesic distance R , and consider the value of volume 
Ul x (R) which gets within the ball d[x,y ) < R. The ho¬ 
mogeneity of the foam means that < Ul x (R) >= D(I?), 
i.e., does not depends on the position of the starting point 
for geodesic lines. Then the scaling f l(R) oc R Dh defines 
the Hausdorff dimension of space Dh- By the construc¬ 
tion it is clear that in 4d gravity Dh < 4. Indeed, the 
starting point x together with the omnidirectional jet of 
geodesic lines define the extrapolating reference system 
(exactly as we use in astrophysics) for which the coor¬ 
dinate volume scales with D = 4. The same behavior 
works for simple topology spaces. However, in the pres¬ 
ence of wormholes for sufficiently big distances R some 


of geodesics return and start to cover the same (already 
covered) physical region Ul(R). Therefore, for sufficiently 
complex topologies (foam) the dimension is always less 
than 4 (i.e., we may state that Dh < 4). 

It is remarkable that we see on the sky exactly such 
a picture. Indeed, the light is too scattered upon prop¬ 
agating through the wormhole throats and it forms the 
diffused background @]. Therefore, in the homogeneous 
Universe galaxies are good tracers for the actual (physi¬ 
cal) volume which exactly demonstrate the fractal behav¬ 
ior N ( R ) oc R d with D ps 2 up to distances R ~ 200Mpc 
[S]. In particular, this reflects the long standing puz¬ 
zle of missing baryons. The nontrivial complex topology 
leads to the modification of Newton’s law [g] which is 
interpreted as the dark matter phenomenon and which 
perfectly fits the observations [7f. 

Thus, we see that lattice quantum gravity suggests us a 
new cosmological scenario. The inflationary stage in the 
past has enormously stretched all scales and the fractal 
quantum spacetime foam structure has been tempered. 
We may say that together with metric perturbations gen¬ 
erated from quantum fluctuations the inflation generates 
topological perturbations from quantum spacetime foam 
as well. Such a structure represents a homogeneous and 
isotropic space filled with a gas of wormholes which repre¬ 
sents the initial conditions for the standard cosmological 
model. In other words, such a space is homogeneous, 
while the topological structure is not (its mean value is 
homogeneous and isotropic but possesses fluctuations). 

Expecting the question on the stability of such worm- 
holes we point out to the well known fact that the non¬ 
trivial topological structure (the gas of wormholes) is con¬ 
sistent with the homogeneity of space. Indeed, such a 
structure can be realized as space of a constant negative 
curvature (by means of the cut and paste technique in 
the Lobachevsky space). In the three dimensions this 
however assumes that the simplest wormhole has throat 
sections in the form of a torus or a sphere with a single 
handle, since the sphere does not possess the metric of 
a constant negative curvature (at least one point will in¬ 
clude a singularity). If the handle on such a sphere is suf¬ 
ficiently small, it will look almost as a sphere (e.g., as we 
know the simplest model of a horse is given by the spher- 
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ically symmetric horse in vacuum; In this respect all our 
attempts to work with spherically symmetric wormholes 
look like this and they may have in the first place the 
methodological interest). Nevertheless, upon averaging 
out over orientations of the torus the spherical symme¬ 
try restores and we may expect that spherical wormholes 
reproduce correctly some basic features. We present the 
simplest exact example of the constructing of a stable 
wormhole in the appendix. 

In general, spacetime foam assumes that entrances into 
wormholes may be separated by time-like and space-like 
distances equally. However it is clear that due to rapid 
inflationary process the time-like separation does not sur¬ 
vive (it remains of the order of Ipi) while space-like dis¬ 
tances have enormously change oc e HAt , where At is the 
duration of the inflationary stage. According to [3] the 
volume of such a space scales as V(R) oc R dH with the 
Hausdorff dimension between 1 < dp < 3. At very small 
and very large scales (maybe even outside the Hubble ra¬ 
dius) dp = 3, while on intermediate scales the dimension 
changes. 

The difference dp - 3 / 0 we observe now as the 
dark matter phenomenon. Indeed, the simplest estimate 
gives for the Newton’s law F oc where S(R) is 

the surface of the sphere of the radius R which scales as 
S(R) cc R dn ~ 1 . We point out that the value dp — 2 
is in agreement with observations in the range 5Kpc < 
R < 200Mpc. The value dp — 1 gives too strong gravita¬ 
tional coupling and such a topological structures should 
decay (or even be suppressed during the inflation). In¬ 
deed, as we know metric perturbations do not grow dur¬ 
ing the radiation dominated stage. However such a be¬ 
havior works only for Freedman spaces (flat space, sphere 
or the Lobachevsky space). In the presence of a nontriv¬ 
ial topological structures this, in general, is not true, e.g., 
see the first investigation in Q which has reviled an es¬ 
sential (scale-dependent) modification of the equations 
for perturbations. 

The negative curvature of space looks to be forbid¬ 
den by the Doppler picks in AT/T spectrum. However 
the upper bounds on the curvature are model dependent 
(they are based on simple topology spaces). Moreover, 
the apparent value of the curvature is somewhat reduces 
by the ratio V p h(R)/V coor (R), where V p h{R ) ~ R dn is the 
actual or physical volume of space and V coor (R) ~ R 3 is 
the extrapolating coordinate volume. In this estimates 
we should take R = Rh (R-h is the Hubble radius). The 
same ratio describes also the portion of missing baryons. 
In other words, these problems require the further and 
more careful investigation. 

Thus, we see that our Universe can be rather far from 
the simple picture we use. It remains to be isotropic 
and homogeneous but may possess a rather complex local 
topological structure. Our phenomenological description 
based on the standard A CDM model works well enough, 
but meets some small inconsistencies which permanently 
enforce us to add some exotic matter fields (absent in lab 
experiments) or consider different modifications of grav¬ 


ity. However it is rather clear that the polarization of 
matter fields on the fractal topology (bias or topological 
susceptibility) is rich enough and it is capable of explain¬ 
ing all exotic properties observed. 

II. HOMOGENEOUS GAS OF WORMHOLES 

Consider now the simplest model suggested by us ear¬ 
lier j| to demonstrate that the topological bias allows 
to mimic dark matter phenomena. As it is explained in 
the previous section stable wormholes have throat sec¬ 
tions in the form of torus (e.g., see Fig.3). Therefore 
the model based on spherical wormholes (which are not 
stable) has in the first place the methodological inter¬ 
est. Nevertheless, it contains all basic general qualita¬ 
tive features (topological permeability). Moreover, upon 
averaging over orientations of tori (wormhole throat sec¬ 
tions) the spherical symmetry of wormholes restores and 
we may expect that even some quantitative features will 
be correct. 

When the gravitational field is rather weak, the Ein¬ 
stein equations for perturbations reduce to the standard 
Newton’s law 

= 47tG ^5p + > 

here a is the scale factor of the Universe, Sp and Sp are the 
mass density and pressure perturbations respectively, G 
is the gravitational constant, and A = V 2 is the Laplace 
operator. Therefore, the behavior of perturbations is de¬ 
termined by the Green function 

A G(x, x') = 47 rS(r — r'). 

In the simple flat space the Green function is well known 
Go = — 1/r (or for Fourier transforms Go = —An/k 2 ). In 
the presence of wormholes due to polarization on throats 
the true Green function obeys formally to the same equa¬ 
tion but with biased source (which is the topological bias 
or susceptibility) 

A G{x, x') = 47t (6(r — r') + b{r — r')). 

In the case of a homogeneous gas of wormholes such a 
bias was evaluated first in j6] (see also more general con¬ 
sideration and details in |l0|) and is given by 

b{k) = 2nR^ (u(^) - ^(0)), 

where n is the density of wormhole throats, R is the 
mean value of the throat radius, and z/(fc) is the Fourier 
transform for the distribution over the distances be¬ 
tween wormhole mouths. It is defined as v{X) = 
-= J F(X±, R)RdR, where F is the number density of 
wormholes in the configuration space (due to homogene¬ 
ity X = X + — X _, and X± are positions of wormhole 
entrances in space). This function is normalized so that 
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f is(X)d 3 X = 1. We point out that in a more general 
case (non-spherical wormholes) the bias will have anal¬ 
ogous structure with an appropriate redefinition of the 
dimensional parameter 2nR —»■ l/£ 2 , though some details 
may change. Thus the true Green function will include 
some correction and look like 


G(k) 


— 4:7T 

k 2 (1 - b{k)) 


— 47T 


1+2n n± 


k 2 



In order to get the logarithmic correction (observed in 
galaxies) the distribution v{k) should give the behavior 
b(k ) ~ k a with a ~ 1. 

Consider now particular examples of the above destri- 
bution. Let all distances between wormhole entrances are 
of the same value r 0 , then i/(X) = (47rrg) <5(|X| — ro) 

and we find v{k) = j v{X)e~ ikx d 3 X = Ofg*, which 
gives 


_ A.'jr 

b(k) = —2nR— ( 1 - 


sin kro \ 
kr 0 ) ' 


For small kro -C 1 we find 

K k ) = Y nRr o + | (krof + ...^ . 


actual distribution of wormholes and restore the true 
Green function 

G ==^(C + Bk a + ...). 

For observational needs as an empirical Green function 
(or generalized susceptibility) we may suggest the expres¬ 
sion 

_ -47 r 

emp k 2 (1 + (fcr 0 )-“) 

which at small scales (kro 1) gives the standard New¬ 
ton’s law, while at large scales kro -C 1 transforms to the 
fractal law. We also point out that the logarithmic cor¬ 
rection (observed in galaxies) corresponds to the value 
a ~ 1 and ro ~ 5 Kpc. 


III. APPENDIX 

Consider now the cut and paste technique on the 
Lobachevsky space which allows to construct stable 
wormholes. Consider the upper complex half-plane, 
see Fig.l, which gives the model of 2D Lobachevsky 
space. The metric on the half-plane has the form dl 2 = 


The first constant simply renormalizes the gravitational 
constant, while next terms define corrections to the 
NewtOn’s law. At first look such a decomposition rep¬ 
resents a rather general situation. Indeed, Consider an 
additional distribution over the parameter r 0 with any 
probability density p(x) (/ 0 °° p(x)dx = 1). Then the 
same decomposition works for mean values < sl1 ^/' 0 >= 
]T) ( 2 n+i)! < (kro) 2n > which means that the above ex¬ 
pression works as well with the replacement Tq™ —>■< 
r 2n >. However this is possible only for normal (Gaus¬ 
sian) distributions with stable (finite) momenta. In the 
case of fractal picture this is not true. Indeed, Consider 
a particular fractal distribution with infinite dispersion 
of the type 

v(k) = exp (—A(ik) a — B(—ik ) a ). 
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FIG. 1: The upper complex half-plane. g± are parts of 
geodesic lines which end at the absolute, a is the equal parts 
of the two geodesics, which upon gluing along g± transform 
to the closed geodesic lines as shown on Fig. 2. 


In the case a <2 the dispersion is divergent 
d 2 

cr = — Tj-(lnv(fe))|fc=o c c k a 2 -» oo 

and the above decomposition does not work. However all 
corrections to the Newton’s law can be found from the 
decomposition of the characteristic function v(k) itself 

v{k) - 1 = k a ^2c n k an . 

All the coefficients in the above decomposition reflect de¬ 
viations from the standard Newton’s law and may be in¬ 
terpreted as the presence of dark matter. In the first 
place they should be fixed from observations (the ob¬ 
served distribution of DM) which allow us to define the 


p- ( dx 2 + dy 2 ) . The absolute (infinity) is the axis Ox 
(y = 0 or y —> oo). Geodesic lines are half circles with 
centers on the absolute, or perpendicular to the absolute 
rays. Making the cut along pies-wise geodesic lines as 
shown on Fig.l (solid lines) and gluing along identical 
geodesics g± we obtain the space of a constant negative 
curvature with two identical closed geodesic boundary 
lines a as shown on Fig.2a. Again gluing along closed 
geodesic circles a we get the Lobachevsky space with a 
handle on it, as shown on Fig.2b. It is important that 
such a space is the space of the constant curvature. It is 
also clear that repeating such a procedure we may insert 
an arbitrary number of handles. 

Now using the axial symmetry of Lobachevsky space 
we add to the above 2d wormhole an angle 0 < <fi < 2n, 
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i.e., rotating the above construction around the axis 
Oz as shown on Fig.3, we get the simplest stable 3d 
wormhole. We repeat again that all these construc- 



FIG. 3: The ball is the 3D Lobachavsky space. Rotation of 
its section around the axis Oz transforms the geodesic circles 
a on Fig.2a into the tori. Insides of tori are removed. Upon 
gluing by the surfaces of the tori we get the simplest stable 
wormhole (handle) in 3D Lobachevsky space. 


FIG. 2: a) Upon gluing along g± we get two identical closed 
geodesics a. Internal region of the circles a is removed, b) 
Upon gluing along a we get the handle in the Lobachevsky 
plane. Geodesics g± transform into a single continuous 
geodesic line. 


tions represent spaces with a constant negative curva¬ 
ture. Their subsequent cosmological evolution is gov¬ 
erned by the Freedman equations (i.e., the metric takes 
the form ds 1 2 3 4 5 = dt 2 — a 2 {t)dl 2 , where dl 2 corresponds to 
the Lobachevsky space with a set of wormholes described 
above). 
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